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ABSTRACT 

We compare the properties of galaxy groups extracted from the Updated Zwicky 
Catalogue (UZC) with those of groups extracted from TV-body simulations of the local 
Universe, in a ACDM and a rCDM cosmology. In the simulations, the initial conditions 
of the dark matter density field are set to reproduce the present time distribution of 
the galaxies within 80/i~^ Mpc from the Milky Way. These initial conditions minimize 
the uncertainty originated by cosmic variance, which has affected previous analyses 
of this small volume of the Universe. The simulations also model the evolution of the 
photometric properties of the galaxy population with semi-analytic prescriptions. The 
models yield a galaxy luminosity function sensibly different from that of the UZC and 
are unable to reproduce the distribution of groups and their luminosity content. The 
discrepancy between the model and the UZC reduces substantially, if we redistribute 
the luminosity among the galaxies in the simulation according to the UZC luminosity 
function while preserving the galaxy luminosity rank. The modified ACDM model 
provides the best match to the UZC: the abundances of groups by harmonic radius, 
velocity dispersion, mass and luminosity are consistent with observations. We find 
that this model also reproduces the halo occupation number of groups and clusters. 
However, the large-scale distribution of groups is marginally consistent with the UZC 
and the redshift-space correlation function of galaxies on scales larger than 6h~^ Mpc is 
still more than S-cr smaller than observed. We conclude that reproducing the properties 
of the observed groups certainly requires a more sophisticated treatment of galaxy 
formation, and possibly an improvement of the dark matter model. 

Key words: methods: miscellaneous - galaxies: clusters: general - galaxies: formation 
- cosmology: miscellaneous - dark matter - large-scale structure of Universe. 



1 INTRODUCTION 

The formation and evolution of galaxies is one of the 
major challenges of cosmology. Within the framework 
of the hierarchical structure formation by gravitational 
instability, galaxies properties are strongly affected by 
their environment. Most galaxies reside in groups and 
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that these systems can be more effective than clusters at 
shaping the properties of galaxies, namely morphology, 
luminosit y, color, star fo rmation rate and cold gas content 
(see also iKodama et"al.i 2001. : .Ellingson. .2004 : .Verheiieiil 
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Galaxy groups trace the large-scale distribution of 
galaxies (e.g. lBerlind et al]l2006^ . and thus describe the dis- 
tribution of the optical light from sub-megaparsec scales to 
the largest scale probed. Models of galaxy formation should 
be able to reproduce the properties of galaxy groups. For an 
appropriate comparison of models with observations, how- 
ever, the surveyed volume of the real groups must be suffi- 
ciently large to be representative of the Universe. Over the 
last years, the compilation of deep galaxy redshift surveys 
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has provided galaxy group catalogues with an increasing 
number of systems (iRamella et al.l[l999l : lEke et all l2004al : 
iBerlind et al.l2o'o^ and references therein) . These catalogues 
are deep enough to average out the large-scale structure fluc- 
tuations and obtain a robust estimate of the light distribu- 
tion. 

In current models, large-scale structure forms by the 
gravitational collapse of dark matter, which dominates 
the dynamics of galaxy systems. Comparing models with 
data therefore requires a valuable treatment of the dynam- 
ics of the baryonic matter within the dark matter halos. 
IPiaferio et alJ l)l99!j) were the first to attempt to compare 
real redshift surveys with mock surveys where galaxies were 
formed and evolved with semi-analytic prescriptions applied 
to merger trees of dark matter halos extracted from A^-body 
simulations. They compared the model with the northern 
secto r of the Updated Zwicky Catalog (UZC) redshift sur- 
vey ifPalco et alJll999fl . This comparison resulted in a gen- 
eral agreement between the now standard Cold Dark Matter 
model with a non-zero cosmological constant (ACDM) and 
observations, although the mean luminosity of the simulated 
groups was smaller than observed. A further topological 
analysis of the mock redshift surveys and the UZC quantified 
how the large-scale structure was much fuzzier in the moc k 
catalogues than in the UZC ijSchmalzing fc Diaferiol [200(1) . 
Both disagreements were tentatively attributed to cosmic 
variance, responsible for the lack, in the mock surveys, of a 
coherent structure as thin and wide as the Great Wall. The 
UZC is only 150/i~^ Mpc deep and indeed its volume is not 
large enough to suppress the fluctuations of the large-scale 
structure. Therefore a fair comparison of the model with 
the UZC requires a not-random choice of the volume in the 
sim ulations. 

iMathis et al.l (|20o3) perform A'^-body simulations where 
the initial conditions are constrained to reproduce the large- 
scale di stribution of galax ies observed in the IRAS 1.2 Jy 
survey jPisher et al.lll99,'Tl) within 80h~^ Mpc of the Milky 
Way. In their simulations, iMathis et alJ ll2002f) also include 
semi-analytic recipes to form and evolve galaxies. Thus, 
these simulations seem appropriate to test whether the cur- 
rent galaxy formation models yield the observed clustering 
and luminosity properties of groups once the large scale dis- 
tribution of galaxies faithfully mirror the real one. This pa- 
per is devoted to this test. 

Sect. 2 summarizes the properties of the constrained 
simulations; sects. 3 and 4 describe our real and mock galaxy 
redshift surveys and our group catalogues in redshift space. 
In sect. 5 we investigate the large-scale distribution of groups 
and in sect. 6 we discuss the distribution of light among and 
within groups. We conclude in sect. 7. 



2 THE SIMULATION OF THE LOCAL 
UNIVERSE 

The constrained simulations we use here combine the A'^- 
body technique with the semi-analytic approach for the for- 
mation of galaxies to predict the evolution of both the pho- 
tometric properties of galaxies and their p hase space coor- 
dinate s, fo llowing a strategy pionereed by iRoukema et alJ 
l|l997l) and lKauffmann et alJ lll99gll. 

The simulations were run by IMathis et alJ ll2002f) who 



provide extensive details in their original paper. They in- 
vestigate two variants of a flat cold dark matter (CDM) 
universe: a ACDM model, with cosmological density param- 
eter fio = 0.3, cosmological constant A — 0.7, and Hubble 
constant Ho — 70 km s~^Mpc~^, and a rCDM model, with 
fio = 1 and Ho = 50 km s~^Mpc~^. For both models, the 
normalization of the power spectrum of the initial density 
perturbations is imposed by the abundance of galaxy clus- 
ters at the present time: erg = 0.9 and 0.6 for the ACDM 
and rCDM models, respectively. 

The initial conditions at z = 50 are set to reproduce 
the observed present-da y galaxy density fi eld, provided by 
the IRAS 1.2 Jy survey (iFisher et alJll994l: 1995), within a 
sphere of radius 80 Mpc centered on the Milky Way 
and smoothed with a Gaussian filter with one-dimensional 
dispersion 5 Mpc. The constrained initial displacement 
field for the dark matter particles for all wavenumbers in 
the range fco = 2n/L and 64fco, on a simulation cube of side 
L = 240/ t~^ Mpc, is generated w ith the Hoffman-Ribak al- 
gorit hm jHoffman fc Ribak|[l99ll 1992; ICanon fc HoffmanI 
Il993r . The high-frequency field in the wavenumber range 
64fco — 343fco is unconstrained. Therefore, the distribution 
of matter on scales smaller than 5 Mpc is unrelated 
to the real Universe. The simulation contains two popula- 
tions of particles: the high-resolution particles, which will 
end up within 80 Mpc from the center of the simulation 
box (the location of the Milky Way) by 2 = 0, and more 
massive low-resolution particles, which will remain outside 
this region. The simulation box contains ~ 5 x 10^ high- 
resolution particles with individual mass 0.36 X 10^°/i"^Mo 
and 1.2 x lO^°/i"^M0 in the ACDM and rCDM model, re- 
spectively. 

The equations of motion of the dark matter par- 
ticles are int egrated with the tree-code GADGET 
jSpringeL Yosh ida & White ,2001,). The halo merger trees 
are then extracted from the simulation and provided as 
input to the semi-analytic code used to form and evolve 
galaxies wi thin the dark matter halos. Apart from a few 
differences, 'Mathis et aU l)2002h use the same recipes that 
Kauffmann et al. ( 1999) applied to the GIF simulations. The 
simulations have enough mass and spatial resolution to fol- 
low the formation of all galaxies brighter than the Large 
Magellanic Cloud that will end up, hy z — 0, within the 
nearly spherical high-resolution region of radius 80 Mpc 
centered on the Milky Way. The luminosity resolution of the 
simulations are Mb = — 15.48 -I- 5 log h and — 16.94 -I- 5 log h 
for the ACDM and r CDM models, respec tively. Note that 
the limits mentioned in lMathis et al.l l)2002h differ from those 
provided on their web site. Here we adopt the latter. To these 
limits, the ACDM and rCDM models contain 189122 and 
194243 galaxies in the high-resolution volume. 



3 MOCK AND REAL REDSHIFT SURVEYS 

We extract mock redshift surveys from the si mulations to 
perform a detailed comparison with the UZC ijFalco et alJ 
..1933) ■ This catalogue is currently the magnitude-limited 
redshift survey of the local universe with the widest sky 
coverage in the optical ba nd. The UZC has b een exten- 
sively inv estigated :1Sc hmalzing fc Diaferiol ll200 d) anal vze its 
topology. IPadmanabhan. Teemark fc HamiltonI il200ll) com- 
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Figure 1. The B-band luminosity functions in the simulations. The vertical lines show the luminosity resolu tion. The dotted lines are 
the fi ts to the Schechter function {Table0. The solid lines are the Mz„-band luminosity function of the UZC iMarzke. Huchra fc Gelle3 
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performed a detailed comparison of the GIF simulations 
iKauffmann et al. 1999) with the CfA2N, the subsample of 
the UZC galaxies in the northern galactic emisphere. 

Here, we consider the region limited by —2.5° ^ <5i950 ^ 
50° and 8'' ^ aigso sj 17'' in the North Galactic Cap (here- 
after NGC) and by 20'' ^ aigso < 4*^ in the South Galactic 
Cap (hereafter SGC). This survey region is 98 per cent com- 
plete down to the Zwicky magnitude mzw = 15.5 (Falco et 
al. 1999). We also apply a Galactic cut —13° ^5 6 ^ 13° and 
further discard the southern Galactic region with Q1950 ^ 3'' 
to reduce extinction problems (Padmanabhan et al. 2001; 
Ramella et al. 2002) . Finally, we only consider galaxies with 
redshift cz < 8000 km s~^, because the simulated high- 
resolution region is limited to 80h~^ Mpc from the Milky 
Way. The redshift surveys we extract cover roughly a quar- 
ter of the sky and contain more than 7000 galaxies. 

To extract mock redshift surveys from the simulation 
box, we place an observer at the center of the box where 
the Milky Way is located and assign celestial coordinates to 
the simulated galaxies. Each galaxy in the simulations has 
radial velocity cz — v ■ r/\\r\\ + \\r\\, where v is the galaxy 
peculiar velocity and r is the galaxy position vector, in units 
of km s^^ . Further details on the construction of the samples 
are given below. 



Table 1. Luminosity function parameters. For the UZC A/* is in 
the Zwicky system. 





a 


- 5 log h 


ij>*/h^Mpc ^mag~l 


UZC 


-1.00 


-18.80 


4.00 X 10-2 


ACDM 


-1.82 


-21.08 


8.39 X 10""' 


rCDM 


-1.08 


-20.29 


1.28 X 10-2 



3.1 The SALF catalogues 

3.1.1 Mock catalogues 

The UZC is magnitude limited to the Zwicky magni- 
tude mzw = 15.5, whereas the mock catalogues con- 
tain galaxies with B-band apparent magnitudes ma = 
Mb + 25 + 51og(r//i-^Mpc). Therefore, we need to trans- 
form the Zwicky magnitude limit into the simulated B- 
band. Detailed calibrations of observed samples seem to 
imply TOzw ~ ru B, but with a rather large 1-cr scat- 
ter of ~ 0.3 mag llHuchral Il976l: ISothun fc Cornel Il990l : 
iMarzke. Huchra fc Gelleil Il994) . Moreover the luminosity 
functions of the galaxies in the simulations and in the 
UZC are sensibly different, as we discuss below, and assum- 
ing naively mzw = would not provi de re liable results. 
Therefore, following IDiaferio et al.l lll999ll and lMathis et all 
i2002h . we adjust the magnitude cut-off of the simulated 
surveys to obtain a number of galaxies in the mock surveys 
comparable to the number in the UZC. 

This constraint imposes the magnitude limits 

(i) ms = mzw + 0.7 = 16.2 in NGC 

(ii) TUB = mzw + 0.4 = 15.9 in SGC 

for the the ACDM model, and 
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(i) ms = mzw - 0.6 = 14.9 in NGC 

(ii) ms = mzw - 0.9 = 14.6 in SGC 

for the rCDM model. Unlike iMathis et alJ ll2002h . we there- 
fore use different transformations for NGC and SGC. These 
different limits are consistent with the 0.3 mag scatter men- 
tioned above and could be due to differences in the photo- 
metric zero-point across the different regions of the sky. A 
few of the many claimed systematic errors in the Zwicky cat- 
alogue are supported by clean observational evidence and a 
unique interpretation of this discrepancy between NGC and 
SGC is still lacking (|Marzke. Huchra fc Gelleil Il994^ . We 
note that the transformations for NGC are similar to those 
adopted by ^iaforio ct al. ( 1999,) in the GIF simulations 
w hich contained semi- analytic prescriptions similar to those 
of lMathis et al.l i2002h . In the UZC and mock catalogues we 
also exclude galaxies with redshift cz < 500 km s~^ and 
mzv < 10 to avoid bright objects close to the Milky Way. 



3.1.2 The UZC 

In order to make a proper comparison with the simulations, 
in the UZC we have to exclude galaxies that are fainter 
than the luminosity resolution of the corresponding simula- 
tion. This is especially important for NGC, where the nearby 
Virgo cluster has many faint galaxies that cannot be resolved 
by the simulations. 

Therefore we compile two UZC catalogues, one where 
we exclude all galaxies fainter than the ACDM luminosity 
resolution (Ms = -15.48 + 51og/i: UZC-A-SALF) and an- 
other where we exclude all galaxies fainter than the rCDM 
luminosity resolution (Afs = —16.94 + 51og/i: UZC-r- 
SALF). 

In principle, the luminosity resolution cut should be 
done in the Zwicky and not in the B-band magnitude. As we 
have discussed above, the transformation between the two 
systems is uncertain and possible differences in the zero- 
points of NGC and SGC make the transformation more 
complicated, because different cuts would be required for 
different regions of the sky. Fortunately, the luminosity res- 
olution cut only concerns the nearby and faintest galaxies 
and a change of ±0.3 mag in the luminosity resolution limit 
reflects into a difference of ±170 galaxies, at most, in the 
final UZC number of galaxies. Furthermore, in the UZC the 
luminosity resolution cut applies in redshift space, whereas 
the luminosity resolution limit of the simulation is in real 
space. When building mock redshift surveys, the effect of 
peculiar velocities blurs such a limit. Given such uncertain- 
ties, it is reasonable to use the simulated B-band limits to 
cut the UZC galaxies in the Zwicky magnitude. 



3.2 The UZCLF catalogues 

The properties of the galaxy distribution in a magnitude- 
limited redshift survey depend both on the distribution of 
galaxies in real space and on their luminosity function. Fig- 
ure shows the galaxy luminosity function in the two sim- 
ulations, with their best fits to a Schechter function (Table 
. The model functions are sensibly di fferent from the lum i- 
nosity function of the UZC derived bv lMarzke et alj il994l) . 
In order to separate the role of the large-scale distribution 



of galaxies from their luminosity function in our compari- 
son between the models and the UZC, we replace the lumi- 
nosities of galaxies predicted by the semi-analytic procedure 
with l uminosities extract ed from the UZC luminosity func- 
tion ( Marzke et al.lll994^ . 

For each model we find the new luminosity resolution 
limit iiim which satisfies the relation 

poo 

V / <t)vzc(L)AL = N (1) 

where 4'vzc{L) is the UZC luminosity function (see Ta- 
ble , A'^ is the number of simulated galaxies more lumi- 
nous than the semi-analytical luminosity resolution and V is 
the high-resolution volume of the simulation box. We then 
randomly sample A'^ luminosities in the range [Liiin,±oo] 
which distribute accordingly to (jjvzciL) and assign these 
new luminosities to the simulated galaxies while preserv- 
ing the luminosity rank predicted by the semi-analytic pre- 
scriptions. The new luminosity resolutions Liim are Afzw = 
-15.85 ±5 log /i and Mzw = -15.78 + 5 log /i for the ACDM 
and rCDM models, respectively. The luminosity function of 
the simulated galaxies now perfectly matches <j)uzc{L) and 
the only galaxy property imposed by the model is the galaxy 
luminosity rank. 

Even if the luminosities are now constrained to match 
the UZC luminosity function, a magnitude limit ~ 15.5 
for the mock catalogues does not necessarly return the same 
number of galaxies as observed. In fact, both the afore- 
mentioned uncertainties in the Zwicky magnitudes and the 
large-scale distribution of galaxies in the small UZC volume 
can affect the galaxy number of a magnitude-limited survey. 
Even if the simulations we use are intended to minimize the 
cosmic variance uncertainty by reproducing the gross fea- 
tures of the large-scale structure, the total number A*' of 
the simulated galaxies still depend on the semi-analytic pre- 
scriptions. Therefore, we adopt the following limiting magni- 
tudes for the mock catalogues, in order to recover a number 
of galaxies approximatevely equal to that of the UZC : 

(i) miin, = 15.3 in NGC 
(u) miini = 15.0 in SGC 

for the the ACDM model, and 

(i) miim = 15.5 in NGC 

(ii) mum = 15.1 in SGC 

for the rCDM model. As for the SALE catalogues, the differ- 
ence between NGC and SGC is ~ 0.3 mag. All the apparent 
magnitudes of the galaxies are then increased by the differ- 
ence between the appropriate and 15.5 in order to have 
the same zero-point of the photometric system in NGC and 
SGC. We now have mock catalogues with the same number 
of galaxies and the same luminosity function as the UZC, 
but with the large-scale structure and the velocity field pre- 
dicted by the underlying cosmological model. 

As for the SALE catalogues, we also make different cat- 
alogues from the real UZC, where we exclude all galaxies 
fainter than the appropriate luminosity resolutions: we thus 
have a UZC-A-UZCLE and a UZC-r-UZCLE catalogue. 

Table 121 lists the number of galaxies in each catalogue. 
Figures 121 show the UZC and the mock redshift surveys. 
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Figure 3. The distribution of tlie galaxies (dots) in the magnitude-limited UZC-A-UZCLF and in the mock ACDM-UZCLF redshift 
surveys. 



4 THE GROUP IDENTIFICATION 

We identify groups in redshift space in tlie UZC and in 
our mock catalogues witti the friends-of-friend s (FOF ) al- 
gorithm described in iRamella. Pisani fc Gelleil (119971) . We 
use the linking parameter Vo — 350 km s~^ and the num- 
ber density contrast 5n/n — 80 at the fiducial velocity 
Vf = 1000 km s~^. These parameters minimize the fraction 



of interlopers and provide the best est imates of the proper - 
ties of groups identified in real space ijPiaferio et aljll99tj l. 
Nevertheless, these linking parameters still yield group cat- 
alogues where ~ 40 per cent of the triplets and 20 per 
cent of the groups with four or more members are acciden- 
tal super position s of galaxies (spurio us groups). Therefore , 
following iMahdavi et al.1 boodt) and iRamella et all J2002l) . 
we restrict our catalogues to groups with N ^ 5 members 
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Figure 5. The distribution of tlie galaxies (dots) in tlie magnitude-limited UZC-t-UZCLF and in the mock rCDM-UZCLF redshift 
surveys. 



when we derive the physical properties of groups. However, 
to limit the shot-noise, we include all groups with N 3 
members when we compute the two-point correlation func- 
tion of groups. 

We consider groups with mean velocities in the range 
500 < cz < 7000 km s~^. The lower limit avoids groups 
close to the Milky Way and the upper limit excludes groups 
which are too close to the edge of the mock survey. Figures 



EJini show the location of the groups in the UZC and in the 
model surveys. 



Galaxy Groups in the Local Universe 7 





-2.5= 50 O" -2.5° SiS 50.0'' 

118 groups 63 groups 



Figure 6. The distribution of groups (diamonds) in the UZC-A-SALF and in the mock ACDM-SALF redshift surveys. 
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Figure 7. The distribution of groups (diamonds) in the UZC-A-UZCLF and in the mock ACDM-UZCLF redshift surveys. 



5 GROUPS AS TRACERS OF THE 
LARGE-SCALE STRUCTURE 



Figures H - |K1 show that both the ACDM and the rCDM 
models reproduce the gross features of the local Universe, as 
imposed by the initial conditions. However, neither model 
has structures as sharply defined as in the real Universe. 
In particular, the rCDM model has much looser structures 



than the ACDM model and the UZC, and this reflects into 
much lower fractions of galaxies in groups (see Table I^J. 

ISchmalzing fc Diaferid i200Cl) quantified this disagree- 
ment in mock redshift surveys extracted from the GIF sim- 
ulations. They used the Minkowski functionals to show that 
the UZC has a substantially larger degree of planarity than 
the models. They concluded that cosmic variance can be 
responsible for the inability of the models at reproducing 
large two-dimensional structures like the Great Wall, which 
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Figure 8. The distribution of groups (diamonds) in the UZC-r-SALF and in the mock rCDM-SALF redshift surveys. 
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Figure 9. The distribution of groups (diamonds) in the UZC-r-UZCLF and in the mock rCDM-UZCLF redshift surveys. 



is the main large-scale structure feature in the UZC. Unfor- 
tunately, the Great Wall, which is at ~ 7000 km from 

the Milky Way, is just on the border of our mock catalogues 
extracted from the constrained simulations. The mock cat- 
alogues do indeed show a concentration of galaxies at the 
location of the Great Wall, but this concentration appears 
substantially looser. 

We compute below the two-point correlation functions 
of groups. However, this function is not exhaustive of the 



large-scale structure description when the density field is not 
Gaussian, as in the case of the small UZC volume. Therefore, 

wo provide here an alternative simple test of clustering. 

Consider the space S — {a,cz), the two-dimensional 
projection of the redshift space {a,5,cz). One can consider 
the galaxies as a random sample of an underlying density 
distribution N2D in this space. If galaxy groups are fair trac- 
ers of the large-scale distribution of galaxies, we expect a 
positive Kolmogorov-Smirnov (KS) test between the distri- 
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Table 2. Total number of galaxies in the real and mock redshift 
surveys. 



SALF UZCLF 





UZC 


ACDM 


UZC 


ACDM 


NGC 


5114 


5037 


5006 


4979 


SGC 


2698 


2667 


2692 


2566 




UZC 


rCDM 


UZC 


tCDM 


NGC 


4558 


4621 


5032 


5026 


SGC 


2623 


2699 


2693 


2691 



Table 5. Comparison between the large scale distribution of the 
UZC and the simulated groups. 



SALF UZCLF 



ACDM rCDM ACDM rCDM 
NGC 3e-4/0.06 9e-8/0.03 2e-5/0.01 le-4/0.01 
SGC 2e-3/2e-3 3e-4/le-3 4e-3/0.01 0.04/0.01 



Bidimensional KS test in the a - cz plane between the dis- 
tributions of groups in the UZC and the models in the range 
cz G [500, 7000] / [500, 6000] km s"!. 



Table 3. Number of groups and fraction of galaxies in groups. 





SALF 


UZCLF 




UZC 


ACDM 


UZC 


ACDM 


NGC 


114/0.27 


118/0.35 


114/0.26 


85/0.24 


SGC 


69/0.28 


63/0.34 


69/0.28 


51/0.18 




UZC 


rCDM 


UZC 


rCDM 


NGC 


101/0.23 


70/0.16 


114/0.27 


80/0.13 


SGC 


68/0.28 


34/0.11 


69/0.28 


49/0.15 



butions of galaxies and groups in S. Table|l]shows that this 
is indeed the case in the UZC, at a significance level larger 
than 17%. A large significance level is also provided by the 
ACDM model, whereas the rCDM groups trace poorly the 
large-scale distribution of galaxies in most cases. For both 
models the situation moderately improves when the UZC lu- 
minosity function is used, except for the NGC of the ACDM 
model. This disagreement is not due to the difficulty of re- 
producing the Great Wall, because cutting the catalogues at 
6000 km s~^ does not substantially improve the agreement. 

Now, if in the models is comparable to N2D in the 
UZC, we also expect a positive KS test between the models 
and the observations. The KS test applied to the galaxy 
N2D fails because the initial conditions of the simulation 
are not set to reproduce the observed galaxy distribution on 
scales smaller than 5h~^ Mpc. However, groups trace the 
large scale distribution of galaxies fairly enough, and the 
comparison between the N2D of the groups in the UZC and 
in the models is appropriate. Table|K|shows that the KS tests 
fail when we consider the entire survey region; however, the 
significance levels increase when we limit the region to 6000 
km s~^, therefore excluding the Great Wall. 



Table 4. Groups of galaxies as tracers of the large-scale structure. 





SALF 


UZCLF 




UZC 


ACDM 


UZC 


ACDM 


NGC 


0.17/0.65 


0.11/0.07 


0.25/0.54 


0.03/0.08 


SGC 


0.22/0.18 


0.32/0.13 


0.23/0.19 


0.37/0.20 




UZC 


rCDM 


UZC 


tCDM 


NGC 


0.19/0.45 


0.002/0.01 


0.23/0.56 


0.05/0.04 


SGC 


0.29/0.23 


0.02/0.03 


0.23/0.18 


0.21/0.24 



Significance levels of the bidimensional KS test in the a — cz plane 
between the distributions of galaxies and groups in the range cz S 
[500,7000] / [500,6000] kms"!. 



Table 6. Number of galaxies brighter than —19.02 + 51og/i in 
the volume-limited catalogues. 



SALF UZCLF 





UZC 


ACDM 


UZC 


ACDM 


NGC 


1788 


1936 


1788 


1431 


SGC 


941 


1011 


941 


479 




UZC 


tCDM 


UZC 


tCDM 


NGC 


1788 


1823 


1788 


2004 


SGC 


941 


875 


941 


655 



We finally consider groups in volume-limited catalogues 
(Table inj. The ACDM model agrees with the UZC better 
than the rCDM model, which yields extremely low numbers 
of groups (Table 0. 

All these results suggest that the structures in the 
ACDM model are better defined than those in the tCDM 
model, although the agreement with the UZC is not yet sat- 
isfactory. The use of the UZC luminosity function allievates 
some problems but it is clearly not the only relevant ingre- 
dient. 



5.1 The correlation function of galaxies and 
groups 

In the redshift space of the local Universe, the vector Si — 
cZiVi locates a galaxy with redshift Zi and celestial coor- 
dinates Vi = (Qi,(5i). For small angular separations (smaller 
than 50° in our analysis), the components along the line of 
sight (tt) and projected onto the sky (rp) of the pairwise 
galaxy separation s = Si — Sj are 



Table 7. Number of groups and fraction of galaxies in groups for 
the volume-limited catalogues. 





SALF 


UZCLF 




UZC 


ACDM 


UZC 


ACDM 


NGC 


18/0.08 


15/0.09 


18/0.08 


7/0.06 


SGC 


17/0.15 


15/0.13 


17/0.15 


1/0.01 




UZC 


rCDM 


UZC 


tCDM 


NGC 


18/0.08 


3/0.01 


18/0.08 


5/0.02 


SGC 


17/0.15 


0/0 


17/0.15 


1/0.01 
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where 1 = (si + Sj)/2. The two-dimensional redshift space 
correlation function ^{rp,n) measures the excess probabil- 
ity, compared to a Poisson distribution, that a galaxy pair 
has separation (rp,n). To measure 5(rp,7r), we compute the 
distribution of pair separations in the data and in a Pois- 
son realization of the data with the same radial and angular 
selection criteria 

DD{rp,n)nR 



DR{rp,n) no 



(3) 



where DD and DR are the weighted sums (see below) of the 
data/data and data/random pairs with separation rp and 
TT, and no and are the mean densities of the real and 
random samples, respectively. 

The points in the random sample are radially dis- 
tributed according to the selection function 
rM(z) 

tp(z) = 



rlVi 
J —c 



M(z„,i„) 



(j)(M)dM 



(4) 



where the luminosity function 4'{M) has the parameters, 
appropriate to the model considered, listed in Table 
We consider only galaxies and groups in the range cz — 
[500,7000] km s"\ M(2„i„) is the absolute magnitude 
correponding to mum at the fiducial minimum redshift 
c«min = 500 km s~^; M{z) is the absolute magnitude 
corresponding to mum at any given redshift z. We con- 
strain M{z) to be brighter than the luminosity resolution 
of the corresponding catalogue. We adopt the usual as- 
sumption that the group and galaxy se l ection functions 
coincide (iRamella. Geller fc Huchralll990l: iFrederid Il995bl: 
iTrasarti-Battistoni et alJll997l : iGirardi et alJl2000l) . 

In magnitude-limited samples, the pair sums DD and 
DR are weighted to correct for the rapid decrease of the 
galaxy density with distance. In its simplest form, the weight 
is given by the inverse of the selection function ip{z). A more 
appropriate appro ach, that we apply her e, is the minimum- 
variance estimate (iDavis fc Huchrall982l) . This approach re- 
quires the a-priori knowledge of the volume integral of the 
real-space correlation function ^(r). However, the weights 
are robust against different choices of £,{r); moreover, be- 
cause this choice affects each catalogue in the same way, 
the conclusions of our comparison remain unchanged. We 
adopt a power-law correlation function ^(r) = (r/ro)"^, 
with ro = 5.8/i~^ Mpc and 7 = 1.8 for both galaxies and 
groups. 

Figure \Wl shows the S,{rp,n) maps of our SALF and 
UZCLF galaxy catalogues. These figures clearly show that 
the models have a substantially weaker clustering at large 
Tp and small n. This quantifies the lack of large and coher- 
ent structures. The Finger-of-Gods features (at large tt and 
small Tp) axe instead well reproduced, a result that indicates 
that the virial regions of clusters are correctly represented 
in the models. 

We also compute the two-point redshift-space correla- 
tion function ^(s), where the separation s is 



c 



z? + z] 



- 2ziZ 



J cos 6ij 



(5) 



and 6ij is the angular separation of the two objects (galaxies 
or groups) with redshift Zi and Zj. ^(s) is more appropri- 
ate to quantify the clustering of groups, because the small 
number of objects makes ^(rp,7r) too noisy. To improve the 



Table 8. Fit parameters of the redshift-space correlation function 
of galaxies. 





so 
Mpc 


7 


UZC-A-SALF 
ACDM-SALF 


7.3 ±0.3 
5.5 ±0.1 


1.18 ±0.01 
1.63 ±0.01 


UZC-A-UZCLF 
ACDM-UZCLF 


7.3 ±0.2 
5.4±0.1 


1.19 ±0.01 
1.67 ±0.01 


UZC-T-SALF 
tCDM-SALF 


7.1 ±0.2 
5.1±0.1 


1.19 ±0.01 
1.25 ±0.01 


UZC-T-UZCLF 
rCDM-UZCLF 


7.3 ±0.3 
5.0 ±0.2 


1.18 ±0.01 
1.25 ±0.01 



Table 9. Fit parameters of the redshift-space correlation function 
of groups with 3 or more members. 





so 

h-^ Mpc 


7 


UZC-A-SALF 
ACDM-SALF 


6.9 ± 1.0 
6.0 ± 1.0 


1.37 ±0.03 
1.39 ± 0.04 


UZC-A-UZCLF 
ACDM-UZCLF 


6.9 ± 1.0 
7.5 ±0.8 


1.38 ± 0.03 
1.65 ± 0.03 


UZC-T-SALF 
rCDM-SALF 


6.8 ± 1.0 
14.5 ±2.0 


1.40 ±0.03 
1.15 ± 0.02 


UZC-T-UZCLF 
tCDM-UZCLF 


6.9 ± 1.0 
7.4 ± 1.1 


1.38 ± 0.03 
1.49 ± 0.04 



statistics, we include systems with three or more members 
when computing the group ^{s). 

To estimate the uncertainties we use a bootstrap proce- 
dure with 50 resampled samples in the case of galaxies and 
400 in the case of groups. 

Figures llll and IT^ show ^(s) for galaxies and groups. 
The fits £,{s) = (s/soy in the interval 2 - 10h~^ Mpc for 
each model are summarized in Tables |5] and [ 



Our results confirm the measure of ^(s) bv lMathis et alJ 

(^2002) (our Figure [TTI is equivalent to their Figure 16): the 
tCDM model provides a galaxy ^(s) which is systematically 
lower than the UZC ^(s) by a factor 0.6 to 3. The ACDM 
model provides a much better match on scales smaller than 
5/i~^ Mpc, whereas its ^(s) is up to a factor 2 smaller 
at larger scales. We remind that the initial conditions are 
drawn from the IRAS 1.2 Jy survey rather t han from the 
UZC optically-selected galaxies, and indeed iMathis et alJ 
^2002^ show that the 5(s)'s of the model and of the PSCz 
survey do agree on all scales in the ACDM model and on 
scales larger than 5/i~^ Mpc in the tCDM model. To com- 
pare the models with th e PSCz correlation function, how- 
ever, iMathis et all (l2002h assign a far-infrared luminosity to 
each galaxy by matching the simulated infrared luminosity 
function to the observed IRAS luminosity function, a proce- 
dure similar to the reassignment of the optical luminosities 
we apply to compile the UZCLF catalogues. Therefore, the 
fact that this by-hand luminosity assignments can yield sat- 
isfactory results in some bands but not in others suggests 
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Figure 10. The galaxy two-dimensional redshift-spacc correlation function 5{rp, tt). The bold contour indicates ^(rp, n) = 1. The contour 
levels are separated by logarithmic intervals of 0.2 for i{rp,iT) > 1 and by linear intervals of 0.2 for i{rp,n) < 1. The dotted contours 
show ^(r'p,7r) < 0. 



10.0 F 




2 6 10 2 6 10 

s (h"' Mpc) s (h"' Mpc) 



Figure 11. The galaxy redshift-space correlation function. The error bars are 1-a bootstrap errors. 



that the clustering power of the underlying dark matter dis- 
tribution is probably correct, whereas the model describing 
how galaxies acquire their luminosity in the various bands 
has to be improved. 

Unlike the galaxy ^(s), the group ^(s) depends on the 
adopted optical luminosity function. ACDM-SALF model 



is in better agreement with observations than the ACDM- 
UZCLF, which tends to give a higher degree of clustering 
on smaller scales; in any case both models agree with the 
UZC ^(s) within the errors. The rCDM-UZCLF model also 

agrees with observations, whereas the rCDM-SALF is at 
least 2-cr above the UZC. The higher degree of clustering of 
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10.0 



1.0 : 



0.1 




UZC-SALF 

UZC-UZCLF 

ACDM-SALF 

ACDM-UZCLF 



s (h"' Mpc 



Figure 12. The group redshift-space correlation function. The error bars are 1-cr bootstrap errors. 



the rCDM-SALF model is a consequence of the fact that 
most groups reside in the nearby region around Virgo (see 
Figure (HJ . 

Our results show that the amplitude of the UZC group 
correlation function is smaller than or equ al to the am- 
plitu d e of the galaxy cor relation function jjing fc Zhand 
Il988l: iRamella et ai] 119901) . This lower amplitude is con- 
firmed in the mu ch larger group sam ple extracted from 
the 2dF survey JPadilla et alJ Eqq3). Other work sug 



gests instead that groups in srnall i Girardi et al.l l200ut 
iGiuricin era!] 1200 ll: iPadilla et al ]l2oqJr and large samples 
JZandivarez. Merchan fc PadiUal '2003^ are more clustered 
than galaxies. Most of this disagreement however is likely 
to be due to the fact that many poor groups are spuri- 
ous and decrease the amplitude of ^(s) toward the ^(s) of 
field galaxies. Indeed the normalization so systematically in- 
creases wi th the e xclusion of groups with fewer and fewer 
members JPadilla et al 2004). In fact, we also find that, in 
the UZC, the ^(s) of groups with at least five members has 
So = 8.8 ± 2.1/i"^Mpc and 7 = 1.30 ± 0.05, and is larger 
than that of galaxies. 



6 THE LIGHT DISTRIBUTION IN THE 
NEARBY UNIVERSE 

In this section we investigate the luminosity content of in- 
dividual groups. We study the group luminosity function, 
the mass-to-light ratio, the halo occupation number and the 
light distribution among galaxy members. 



6.1 The luminosity function and abundances of 
groups 

To estimate the space density of redshift-space groups we 
weight each group by l/Fmax, where 
, 3 r 



(!)■ 



1 



3CZ5 



(1 + 90) 



(6) 



is the volume sampled by the group. In fact, cz^/Hq is 
the smallest distance between our group distance cut off 



7000 km s~^ and the maximum distance beyond which the 
fifth brightest galaxy member would be fainter than the ap- 
parent magnitude limit miim; moreover, to is the solid angle 
of the catalogue, c the light speed, and qo the deceleration 
parameter which we suppose to be unknown and set equal 
to 0.5. This choice is irrelevant, because, for the small dis- 
tances considered here, the correct value qo = —0.55 for the 
ACDM model yields a Fmax which differs by less than 3 per 
cent from the Fmax estimated with go — 0.5. 

In real space, we consider a group each dark matter 
halo containing N ^ 5 members brighter than the appropri- 
ate luminosity resolution. We compute the total luminosity 
Ltot of a real-space group as the sum of the luminosities of 
the A'' resolved galaxies within the halo, Lgai — Li, and 
the luminosity of the galaxies fainter than the luminosity 
resolution limit Lum 



iV 



-Jfaint noc 



Lun 



(7) 



where is the adopted galeixy luminosity function. 

In redshift space, we use the same equations with the 
only difference that the luminosity resolution limit Lum is 
replaced by the brightest luminosity between the luminosity 
resolution and the absolute magnitude Mu^ = rnum — 25 — 
51og((cz)/_H'o), where {cz)/Ho is the group distance. 

Not surprisingly. Figure 1131 shows that neither the 
ACDM- nor the rCDM-SALF mock catalogues reproduce 
the UZC group luminosity function. With the UZCLF cat- 
alogues the situation improves, but it is still unsatisfactory, 
because the KS and Wilcoxon-Rank-Sum (WRS) tests in- 
dicate that the mock distributions are not drawn from the 
same population of the UZC at very high significance levels 
(Tables inland int. 

Nevertheless the ACDM-UZCLF appears to yield a 
luminosity function closer to the UZC than the rCDM- 
UZCLF. In fact, the rCDM model fails at providing nu- 
merous enough luminous groups, despite the fact that the 
rCDM model produces a number of bright galaxies larger 
than both the ACDM model and the UZC (Figure [T||. How- 
ever, the small amplitude of the galaxy two-point correla- 
tion function (Figure lllll shows that the galaxies are not 
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Figure 13. Cumulative luminosity function. The error bars are 1-a Poisson errors and are only shown for the UZC for clarity. The error 
bars for the models have similar size. 



enough clustered to yield either a sufRciently large fraction 
of galaxies in groups or a sufficiently large number of groups. 
Consequently, the median luminosities of the rCDM groups 
are lower than in the UZC (Table IT?1jl . 

Figure ITH shows the abundances of groups by harmonic 
radius 



Rh 



TT (cz) 

2lh 



N{N - 1) 



EE 

, i— 1 J— i + 1 



tan(6',j/2) 



(8) 



where is the number of members and 9ij are the pairwise 
angular separations, velocity dispersion 



iV 



1 

— ^(czi - {cz) f 



1/2 



and virial mass 



Mvir = 



G 



(9) 



(10) 



To avoid divergences, we set, according to Ramella (private 
communication), tSLn{9ij /2){cz) /Ho = Vij = 0.03ft~^ Mpc 
when Tij < 0.03/i~^ Mpc. The SALF catalogues perform 
very poorly: the medians of these quantities are sensibly 
different from the UZC medians (Table [TI71 and the KS and 
WRS tests indicate that the mock groups and the UZC 
groups are certainly drawn from different parent popula- 
tions (Tables [TTI and [T^ . The only exceptions are the a and 
Mvir distributions of the rCDM models and the M^ir/Ltot 
distribution of the ACDM model. 

The UZCLF catalogues perform substantially better: 
with the exception of the Ltot and Mvir/itot distributions, 
both the ACDM and the rCDM models provide a rather 
good match to the UZC group abundances, with the ACDM 
model yielding the largest significance levels. 

Figure [T^ also shows the number density of the groups 
in real space (dotted lines). For these groups 



Rh — 



N{N -1) 



E E ^ 

i—l J — z + l 



(11) 



N 



1 1/2 



3(iV_ 1) ^ 

i—l 



(12) 



where Vij are the pairwise separations. 



and Mvir = Qa'^Rh/G. 

In the ACDM model, the number densities of groups 
estimated in redshift space are larger than the real-space 
number densities, whereas in the tCDM model the largest 
discrepancy appears in the number density by harmonic ra- 
dius. This result agrees with tha fact that Rh can be overes- 
timated by a factor as large a s two because of the presence 
of interlopers (.Diaferio et al.lil999'l . In the virial mass esti- 
mation, however, this bias is less relevant than the bias on 
cr, because of the different powers that Rh and a have in the 
virial mass relation: on the scale of clusters, where a is al- 
most unbiased, the recovered mass function is in reasonable 
agreement with the real-space mass function. 

A''-body simulations, w here galaxies were iden- 
tified as density peaks ["Nolthenius fc White! Il987l : 
■Moore. Frenk fc White 1993; FYedcric 199 5a b) or formed 
and e volved with semi-analytic prescriptions lDiaferio et alJ 
^9^, show that the FOF algorithm generally returns 
groups with statistical properties comparable to those of 
the real-space groups; however, interlopers and spurious 
groups can severly affect the estimates of the intrinsic 
properties of individual groups and t his bias is more dra- 
matic for groups with fewer member s (iDiaferio et al.lll999l : 
lEke et al.ll2004bl: iBerhnd et alJl2006l) . Figure [Til shows that 
the discrepancy between the average properties of real- and 
redshift-space groups might be larger than usually assumed. 

Recently, group-finder algorithms alternative to the 
FOF have been proposed and might identify real-space 
groups more accurately than the FOF algorithm ( Eke et ^ 
l2004al : lYang et al.ll2005b() . However, the FOF algorithm has 
been traditionally applied to the UZC and we use it here to 
make our results easily comparable to previous work. More- 
over, we are interested in the comparison between models 
and observations when they are both analyzed with the same 
technique, rather than in the determination of the most re- 
liable group finder. It is indeed a relevant result that the 
UZCLF catalogues agree with the UZC remarkably well, de- 
spite the fact that the recover of the average properties of 
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Table 10. Median values of the group properties. 



SALF 



UZCLF 



Rh 

(T 

log M„ir 
log Ltot 

log(Af„,>/Ltot) 

no 



UZC/ACDM/3D 
0.42 / 0.64 / 0.29 
219 / 284 / 230 
13.45 / 13.81 / 13.29 
10.88 / 11.07 / 10.58 
2.59 / 2.70 / 2.70 
0.27 / 0.27 / - 



UZC/tCDM/3D 
0.42 / 0.51 / 0.28 
222 / 217 / 320 
13.48 / 13.54 / 13.57 
10.92 / 10.56 / 10.65 
2.58 / 2.94 / 2.93 
0.27 / 0.82 / - 



UZC/ACDM/3D 
0.42 / 0.44 / 0.29 
219 / 217 / 230 
13.45 / 13.44 / 13.29 
10.88 / 10.64 / 10.49 
2.60 / 2.85 / 2.79 
0.28 / 0.50 / - 



UZC/rCDM/3D 
0.42 / 0.47 / 0.28 
217 / 250 / 320 
13.45 / 13.59 / 13.57 
10.88 / 10.64 / 10.38 
2.60 / 2.96 / 3.20 
0.28 / 0.64 / - 



Median properties of the galaxy groups with more than five members. Rh, f : 
h~^'bAQ, /i-^Lq and /iMq/Lq, respectively. 



Mvir; Ltot and M^ir/Ltot are in units oi h ^ Mpc, km s 



Table 11. Comparison of the UZC and the mock group abun- 
dances: KS test. 



SALF 



UZCLF 





ACDM 


tCDM 


ACDM 


rCDM 


Rh 


10-9 


10-3 


0.15 


0.14 


a 


10-3 


0.10 


0.27 


0.08 


log Mvir 


10-5 


0.11 


0.34 


0.02 


log Ltot 


10-5 


10-11 


10-7 


10-5 


log{Mvir/Ltot) 


0.15 


10-11 


10-4 


10-9 



Significance levels of the KS test for the null hypothesis that the 
UZC and the mock groups are drawn from the same parent pop- 
ulation. 



Table 12. Comparison of the UZC and the mock group abun- 
dances: WRS test. 



The universal value of the mass-to-light ratio is 
(M/L) = pciit^o/j, where pcrit is the critical density of the 
Universe and j = 0*L*r(a -I- 2) is the luminosity density 
derived from the galaxy luminosity function. In our models, 
the real-space relations approaches a universal value at suf- 
ficiently large masses; this value does not correspond to the 
correct Qo, however, because the virial mass Mvir estimated 
in real space is still ~ 30 — 40% larger then the actual mass 
of the dark matter halo jPiaferio et alJ 119991) . In redshift 
space the M/L — M relations do not clearly show this flat- 
tening, but their rather large percentile ranges include the 
real-space relations. The ACDM model matches the UZC 
results rather well, whereas the rCDM model provides too 
large mass-to-light ratios both because it does not provides 
luminous groups and because the real Universe seems to 
have a low value of flo and, consequently, a low value of 
{M/L}. 



SALF 



UZCLF 





ACDM 


tCDM 


ACDM 


rCDM 


Rh 





10-3 


0.16 


0.04 


a 


10-3 


0.39 


0.41 


0.03 


log Mvir 


10-6 


0.25 


0.29 


0.01 


log Ltot 


10-5 





10-7 


10-6 


log{Mvir/Ltot) 


0.01 





10-6 






Significance levels of the Wilcoxon Rank-Sum (WRS) test for the 
null hypothesis that the UZC and the mock groups are drawn 
from the same parent population. 

the real-space groups with the FOF algorithm is only par- 
tially correct. 



6.2 Mass-to-light ratio vs. mass 

Figure \TE\ shows that the mass-to- light ratio A/vir/itot of 
the UZC groups increases with mass, as other sam ples 
of gr oups have already demonstrated ( Marin oni & Hu dsonI 
l2002l: iLanzoni et alJ 120041: lEke et al.l l2004bL 2006). Semi- 
analytic models of galaxy formation show that the depen- 
dence of Mvir /Ltot on mass naturally arises in hier archi- 
cal models iKauffmann et al .ll999llBenson' et al .120001) . and 
extensions of the halo model use this relation as a sensi- 
tive probe to constrain the efficiency of galaxy formation 
and the normahzation of the power spectrum of the pri- 
mordial density flu ctuations erg ivan den Bosch et alJ2003bl: 
iTinker et a^]|2005^ . 



, ^ ^ laxy prope 

ment jYang et alj l20058i: IShethI I2OO5I : ISkibba et alJ lioo^ : 



6.3 Halo occupation number 

The halo occupation number (HON) is the mean num- 
ber of galaxies brighter than a luminosity thresh- 
old in dark matter halos of a given mass; it quan- 
tifies how efficient galaxy format i on is in halos of 
differ ent masses iKauffman n et all Il997t iBenson et 
200fl| Igyy^^iin^rro e^ ^lj|2(J(,)l|: [y^n (j^^n B9^9h ^l.|l29(^jaf 
iBerlind et al.l2003l : lKravtsov et al.l2004l:IZheng et al.l2005r . 
With appropriate extensions, one can use the HON to 
model the dependence of galax y properties on environ 

^_rig _eLal 
ISheth et alJl200l 

Estimating the HON of real systems requires accu- 
rate photometry and accurate determination of the to- 
tal mass and galaxy membership. To avoid these difficul- 
ties, a common approach is to assume an analytic form 
for the HON and adjust its parameters to mat ch the ob- 
servcd galaxy two-point correlation function (Jing et al 

; Ma gliocchetti & Porciani 2003; Abazaiian ot al. 200:_ 

Zehavi et aliboOfill) . Alternatively, one can directly use the 
galaxy luminosity function to esti mate the expected number 
of galaxies in th e cluster volume iMarinoni fc Hudsonll20d^ : 
iLin et aljliool . 

These approaches have provided an indirect comparison 
of the observed HON with A'^-body or semi-analytic simu- 
lations. A direct comparis on requires galaxy counts for a 
sample of galaxy systems (iKochanek et alJl2003l ILin et alJ 
I2OO4I IPopesso et al.ll200'3) . One of these direct counts was 
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Figure 14. Cumulative number density of groups by harmonic radius Rh, velocity dispersion cr, and virial mass Mvir- The error bars 
are l-cr Poisson errors and are only shown for the UZC for clarity. The error bars for the models have similar size. 




Figure 15. Group mass-to-light ratio vs. virial mass. The Hncs connect the medians of each bin and the error bars, only shown for the 
UZC for clarity, indicate the 0.1 and 0.9 percentiles. The error bars for the models have similar size. The dot-dashed line shows the 
universal mass-to-light ratio in the model. 
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Figure 16. The halo occupation number (HON) in real space. The top (bottom) panels show the ACDM (rCDM) model. In each panel 
lines are for galaxies brighter than Mg — 51og/i = —18.80, and —19.91 (top to bottom) in the ACDM model, and —18.80, and —19.91 
(top to bottom) in the rCDM model. The HON values smaller than one are due to haloes containing no galaxies brigh ter than the 
thres hold. The error bars are l-cr Poisson errors. The points with error bars on both axes show the nine CAIRNS clusters jRines et alj 

Eool. 



performed bv lRines et alj (|2004) for the nine CAIRNS clus- 
ters, thanks to their accurate 2MASS infrared photometry 
and spectroscopic uniform sky coverage. Figure lTHl shows the 
CAIRNS cluster results. 

Before comparing the HON of the UZC groups with the 
model groups, we check what model, if any, can reproduce 
the HON of these real clusters. This comparison is more ro- 
bust than the comparison with groups, because the galaxy 
membership in groups is substantially more uncertain than 
in clusters, and the comparison between models and obser- 
vations will be obscured by the large error bars (see Figure 
IT7I below). 

To compare the CAIRNS HON with the HON in the 
simulations we need t o convert the 2MASS A's-band to the 
-B-band. According to l.Iarretl l)200(t) . B — Ks varies between 
2.86 and 3.97, depending on the galaxy morphology; these 
colours have rms uncertainties in the range 0.30-0 . 80. Th ere- 
fore, the luminosity limit imposed bv lRines" et all (|2004) Ues 
in the B-band range [-18.80, -19.91]. Figure [T^ shows the 
HON of the simulations in real space for both luminosity 
hmits. Because early-type galaxies tend to reside in clusters 
and have redder colours, the fainter limit (dashed line) is 
probably the most appropriate for the comparison. 

The SALF catalogues provide a factor of ~ 3 fewer 



galaxies, whereas only the ACDM-UZCLF model provides a 
good match to the CAIRNS clusters. The slight underesti- 
mate still present in this case might originate from the fact 
that, because of the diffi culties at identifying galaxy mem- 
bers, iRines et al.l l)2004l include all the galaxies projected 
onto the cluster. 

When extending the comparison to the scales of groups, 
uncertainties on galaxy memberships and mass estimates 
become larger. Therefore, because the luminosity Ltot ap- 
pears to be a more robust quantity than the virial mass 
Mvir (Eke ct al. 2004b), we measure the HON as a function 
of both Z/tot and Mvir (Figure lTTji . In general the models un- 
derestimate the UZC A'gai — Mvir relation, regardless of the 
luminosity function adopted. On the contrary, the simulated 
A'gai — itot relations agree well with the UZC, particularly 
the ACDM-UZCLF model. The rCDM models reproduce 
the observed slope and normalization, but their relations 
vanish at large Ltot , because their mock surveys do not con- 
tain enough luminous groups. 



6.4 The light distribution within groups 

In the semi-analytic models analyzed here mergings between 
galaxies are allowed only when one of the two galaxies is 
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Figure 17. The halo occupation number (HON) in redshift space for groups with 5 or more members. The top (bottom) panels show 
the ACDM (rCDM) model. The lines are for galaxies brighter than Mg = —19.02 + 51og/i, the luminosity of a galaxy with apparent 
magnitude rrinj^ = 15.5 at the cut off distance 8000km s~^. The dotted lines are the real space HON at the same Mg = — 19.02 + 5 log /i. 
The HON is shown as a function of group mass Mvir (left panels) and group total luminosity Ltot (right panels). The HON values smaller 
than one are due to haloes containing no galaxies brighter than the chosen threshold. The error bars are 1-cr Poisson uncertainties and 
are only shown for the UZC for clarity. The error bars for the models have similar size. 



the central galaxy of the dark matter halo: satellite-satellite 
merging is not implemented. Moreover gas cooling is al- 
lowed only onto the central galaxy and satellite galaxies stop 
forming stars when their internal reservoir of cold gas is ex- 
hausted. These two processes yield a luminosity difference 
ALi2 between the first and the second brightest galaxies in 
a group that can be substantially larger than observed. 

In Figure ITSl we plot AL12 as a function of the group 
total luminosity. The median AL12 is 3.3 X lO^h-^Le for 
the UZC and 9.3 x IO^/i-^Lq for the ACDM-SALF model. 
The AL12 density distribution of the rCDM-SALF model is 
similar to the UZC (its median is AL12 = 3.7 x lO^/i^^L©), 
but the tail of the high-luminosity groups is missing. To 
quantify the itot — AL12 relation, Table IT^ lists the Spear- 
man rank-order correlation coefficients r and their proba- 
bility P: smaller P's indicate stronger correlations. In the 
ACDM-SALF model, the correlation is orders of magnitude 
stronger than in the UZC. By reassigning the galaxy lumi- 
nosity according to the UZC luminosity function the corre- 
lation weakens to the observed level. 

The semi-analytic recipes produce many more bright 
galaxies than observed in the UZC (Figure 0. This fact can 
boost, in some cases, the luminosity difference AL12, as one 
can see by comparing the two ACDM panels in Figure [T51 



Table 13. Correlation between AL12 and the group total lumi- 
nosity. 



SALF UZCLF 





UZC 


ACDM 


UZC 


ACDM 


r 


0.30 


0.63 


0.29 


0.33 


P 


10-4 


10-19 


10-4 


10-4 




UZC 


tCDM 


UZC 


rCDM 


r 


0.31 


0.35 


0.29 


0.36 


P 


10-4 


10-4 


10-4 


10-4 



Spearman rank correlation coefficient (r) and the significance of its 
deviation from zero (P). A smaller value of P indicates a stronger 
correlation. 

We therefore also consider the ratio L2/L1 of the luminosi- 
ties of the two brightest galaxies. The medians are 0.69 for 
the UZC and 0.52 and 0.63 for the ACDM-SALF and the 
rCDM-SALF models, respectively, and confirm that the sec- 
ond brightest galaxy is on average fainter in the models than 
in the UZC. 

A more quantitative comparison is obtained by consid- 
ering the density distributions of the luminosity ratios. Table 
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Figure 18. Luminosity difference between the first and the second rank galaxies in a group as a function of the group total luminosity. 



Table 14. UZC-model comparison of the luminosity ratio be- 
tween the two brightest galaxies in a group. 





SALF 


UZCLF 




ACDM tCDM 


ACDM tCDM 


KS 


10-5 0.01 


0.78 0.67 


WRS 


10-'' 10-3 


0.41 0.38 



Significance levels of the KS and WRS tests for the null hypote- 
sis that the luminosity ratio L2/-L1 between the second and first 
brightest galaxies in the UZC and in the simulated groups are 
drawn from the same parent population. 



1141 lists the significance levels of the KS and WRS tests for 
the null hypothesis that the UZC and the model L2/L1S 
are drawn from the same parent population. Clearly, the 
ACDM-SALF and the UZC groups have different parent 
populations, whereas the 1% significance level of the rCDM- 
SALF is a consequence of the failure of this model to yield 
groups brighter than Ms — 5 log h ~ —22, as we have clarified 
above (Figure The comparison is instead satisfactory 
when the galaxy luminosities are reassigned according to the 
UZC luminosity function (UZCLF catalogues): both models 
now agree with the UZC to a significance level greater than 
38%. 



7 CONCLUSIONS 

We have compared the properties o f groups and clust ers of 
galaxies extracted from the UZC jFalco et alj Il999^ with 
those of systems extracted from mock redshift surveys. We 
have compiled these mock surveys from A'^-body simula- 
tions constrained to reproduce the large-scale distribu tion 
of galaxies in the nearby Universe jMathis et al.ll2002^ . In 
the simulations, galaxies are formed and evolved according 
to a semi-analytic procedure and we are thus able to con- 
strain both the group clustering and the group luminosity 
content. 

By using simulations with constrained initial condi- 
tions, we minimize the possible role played by cosmic vari- 
ance, and differences between real and mock catalogues 
should mostly originate from the galaxy formation recipes. 
Our approach thus differs from previous attempts of com- 
paring simulations with real group catalogues, w here either 
cosmic variance was an issue llDiaferio et al .119991) or the real 
catalog was largo enough that cosmic variance was naturally 
suppressed (Eke et al. 2004a; Bcrlind et al. 2006). 

The gross large-scale distribution of galaxies, including 
the location of the major nearby clusters, is very similar 
in the mock and real surveys, as mostly imposed by the ini- 
tial conditions. However, the simulated large-scale structures 
are not as sharply defined as in the UZC, confirming earlier 
results of mock surveys extracted from unconstrained sim- 
ulations JSchmalzing fc Diaferiol l2000h . This disagreement 
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amplifies substantially when we consider the properties and 
the large-scale distribution of groups. 

The group-finder algorithm strongly depends on the 
galaxy luminosity function. The numerical recipes adopted 
to form and evolve galaxies in the simulations provide a 
galaxy luminosity function which is considerably different 
from that of the UZC. We test that this difference is the 
major responsible for the model inability at reproducing the 
group properties. We assign new luminosities to the simu- 
lated galaxies according to the luminosity function of the 
UZC, while preserving the luminosity rank predicted by the 
model. We thus have new mock catalogues where galaxy lu- 
minosities distribute according to the UZC luminosity func- 
tion (UZCLF), besides the catalogues where galaxy lumi- 
nosities distribute according to the semi-analytic luminosity 
function (SALF). 

Unlike the groups extracted from the SALF and the 
rCDM-UZCLF mock catalogues, the groups extracted from 
the ACDM-UZCLF catalogue have statistical properties in 
satisfactory agreement with observations: group abundances 
by luminosity, harmonic radius, velocity dispersion and mass 
are generally within 3-a errors, or less, from the UZC group 
abundances. These groups also reproduce the UZC relations 
between the group mass-to-light ratio and mass, and the 
galaxy number and mass, namely the halo occupation num- 
ber. Finally, the ACDM-UZCLF groups, similarly to the 
UZC groups and unlike the SALF groups, show a weak cor- 
relation between the luminosity difference between the two 
brightest galaxies in a group and the total group luminosity. 
This result indicates that the two semi-analytic prescrip- 
tions of allowing merging of satellite galaxies only with the 
central galaxy of the dark matter halo and gas cooling only 
onto this same central galaxy produces a too large luminos- 
ity difference between the two brightest galaxies in a group. 

The success in the statistical properties of groups, ob- 
tained by adopting the observed luminosity function, is not 
shared by the large-scale distribution of groups. Specifi- 
cally, the simulated groups in the North Galactic Cap of 
the ACDM-UZCLF catalogue trace the large-scale distribu- 
tion of galaxies at significance levels as much as seven times 
smaller than those of the UZC. This is a consequence of 
the looser large-scale structures in the models: in fact, the 
redshift-space correlation function of galaxies is more than 
3-cr below the observations on scales larger than 6h~^ Mpc 
(Figure imlMathis et a'Lll2002t) and the number of groups is 
25% smaller than in the UZC (Table gj. 

We therefore con c lude t hat (1) the semi-analytic recipes 
used in lMathis et alJ l)2nn2h have seriuos difficulties at dis- 
tributing the luminosity among galaxies correctly; and (2) 
even the most successful ACDM model does not yield cos- 
mic structures as coherent as observed. Whereas the first 
problem might be solved, in principle, with more sophis- 
ticated semi-analytic m odellings (e.g. ICroton et al.l 120061 : 
IWeinmann et al1l2006lJl the second problem appears to be 
more fundamental, although the disagreement with observa- 
tion is less dramatic: in this case, both different dark matter 
models and an appropriate treatment of the interplay be- 
tween gas dynamics and dark matter might be necessary to 
reconcile the model with group properties. 
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